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Abstract 

The exact solution for the scattering of electromagnetic waves on an 
infinite number of parallel demi-planes has been obtained by J.F. Carlson 
and A.E. Heins in 1947 using the Wiener-Hopf method We analyze 
their solution in the semiclassical limit of small wavelength and find the 
asymptotic behaviour of the reflection and transmission coefficients. The 
results are compared with the ones obtained within the Kirchhoff approx- 
imation. 
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1 Introduction 



For most systems classical mechanics corresponds to the ?i — > limit of quan- 
tum mechanics and the nature of classical motion to a large extent determines 
spectral statistics of quantum problems PP, P]. Nevertheless, for certain mod- 
els this correspondence breaks down because limiting classical Hamiltonian has 
singularities and there is no unique way to continue classical trajectories which 
hit them. Though classical mechanics itself is not complete, quantum mechan- 
ics in many cases smoothes singularities and associates with each of them a 
diffraction coefficient which gives a probability amplitude for the scattering to a 
given channel. The simple example of such models is given by diffractive models 
with small-size impurities where the total wave function at large distances from 
an impurity is a sum of the free wave function plus a reflected field. In two 
dimension in a convenient normalization 

D{6f,9i) is the diffraction coefficient for the scattering of the incident plane 
wave of the direction 9i to a reflected plane wave of the direction 9 f . 

Though in classical mechanics such scatters are negligible, their presence 
perturbs greatly (and in a calculable way) quantum mechanical problems (see 
P] and references therein). 

An important class of diffractive systems consists of plane polygonal billirds 
(see e.g. |1]) with billiard corners playing the role of diffraction centers. The 
main difficulty in such type of diffractive models is the impossibility to represent 
everywhere the scattering field as a sum of a free motion plus small corrections 
as in Eq. which forms the basis of the usual scattering theory. 

As an example we present the exact results for the diffraction on a demi-plan 
derived by A. Sommerfeld in 1896 The total wave function for this problem 
obeying the Helmholtz equation 

(A + fc2)^'(z,a;) = (2) 

in the plan (z, x) with the Dirichlet boundary conditions along the semi-infinite 
screen a; > has the following form (in notations of Fig. ^ 



*(r) = e-''''''°''^'f~''^Fi-V2l^cos^^^) 

- e-'''''^°'^'f+'''>F{~V2I^cos^^^), (3) 
where F[u) is the Fresnel integral 

-i7r/4 roc 

F{u) = \ e'' dt. (4) 
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Figure 1: Scattering on the demi-plan (thick Hne). The initial and final rays are 
indicated by thin lines. The dotted lines represent optical boundaries. 



From the expansion of 5'(r ) at large distance it follows (see that the diffrac- 
tion coefficient for this problem is 

Di9f,e,)^^j^ i^. (5) 

cos — cos — 

The important specificity of this diffraction coefficient is that it formally blows 
up in two directions 

9f^TT±e^ (6) 

called optical boundaries (see Fig.^ which separate regions with different num- 
ber of geometrical optic rays which play the role of classical trajectories. As wave 
fields are continuous, the separation of the exact field ^ in the sum of free mo- 
tion plus small refiected fields is not possible in a vicinity of optical boundaries 
which manifests itself as the divergence of the diffraction coefficient 

In practice one can use diffractive coefficient description for all scattering 
angles except intermediate parabolic regions near optical boundaries where the 
dimensionless arguments of F-functions in Eq. Q are of the order of 1 

y- . Sip 

u — Vkr sin — ^ 1, (7) 

and 5ip is the deviation angle from the optical boundaries ©. 

Difficult problems appear when inside these intermediate regions there are 
new points of singular diffractions which is inevitable e.g. for plane polygonal 
billiards. In the semiclassical limit fc — oo certain cases of multiple diffraction 
on singular corners have been computed within the Kirchhoff approximation in 
The long-range nature of the diffraction on sharp corners leads to a strong 
interaction between different singular points and no general formulas for exact 
diffraction coefficients in these cases are known. 
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In Refs. [7] and [S] the exact solution for the reflection of electromagnetic 
wave on an infinite number of parallel semi-infinite metallic sheets (equivalent to 
the scalar Hclmholtz equation with the Dirichlet ,7 and Neumann |5] boundary 
conditions) has been obtained by the Wiener-Hopf method (see also Exact 
expressions for reflection and transmission coefficients (Eqs. (l3Uf) . H31|l . and H32|l 
below) are quite cumbersome. They are given by infinite products of different 
terms each depending non-trivially on the initial momentum k. The purpose of 
this paper is the investigation of these results in the semiclassical limit k —^ oo. 

The plan of the paper is the following. In Section |21 general properties of 
the reflection from periodic chain of demi-plans are discussed. In Section 13 the 
exact expressions for diffraction and transmission coefficients obtained in j7j are 
given and more tractable expressions for the modulus of these coefficients are 
presented. The most interesting case corresponds to the scattering with small 
incident angle when in the intermediate region (UJ there are many singular 
points. In Section 01 two first terms of the expansion of the elastic reflection 
coefficient are derived. In Section[Slthe first terms of expansion of the refiection 
coefficients for the forward scattering in the power of the incident angle are 
calculated and in Section the same is done for large-angle scattering. In 
Section [3 these calculations are generalized for all powers of the incident angle 
provided the condition Q is fulfilled. The exceptional case when demi-plans 
are perpendicular to the scattering plan is treated in Section|Sl In Section|51the 
exact asymptotics are compared with the ones calculated within the Kirchhoff 
approximation and in Section [TUI the summary of obtained results is given. 



2 Generalities 

The configuration of half- planes considered in Refs. [Z], [HI is represented sche- 
matically at Fig. 121 The equally spaced half-planes are parallel to the z-axis 
and their left corners form a scattering plane inclined by an angle a. The initial 
incident plane wave 4*0 (2, x) comes from the left with angle of incidence 9 with 
respect to z-axis. 

*o(2,a;) = e*'=(^™''*+"'"'"^). (8) 

In the far field due to the periodicity of the problem only finite number of waves 
can propagate. The reflected plane waves have the following form 

V[/(''e/)^^^ 2;) = ^iHz cos e'^+x sin e'^) ^ |-g^ 

where the allowed values of reflected angle 9'^ are determined by the usual 
grating equation (see e.g. [S]) 

kd{cosip — cosip'^) — 27rn. (10) 

Here ip and (f'„ are the incident and reflected angles defined as the angles be- 
tween the scattering plane and the incident (respectively, refiected) plane wave 
direction 

ip^a^9, (11) 
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Figure 2: Scattering on half-planes staggered by an angle a. 9 is the angle of 
incidence of the initial plane wave. Dashed line indicates the scattering plane. 
Dotted line is the contour used in the computation of the current conservation. 



and 

ip'„ ^9'^~a. (12) 

d — a/ sin a is the distance along the scattering plane between corners of the 
half-planes, n is an integer. 

In terms of the dimensionless momentum 

Q - - (13) 

TT 

Eq. H1Q(I determines a finite number of possible reflected directions when (p is 
fixed 

cos(/3„ = cos^- — . (14) 

The allowed values of n are restricted by the inequalities followed from this 
expression 

- Q sin^ I < n < Q cos^ |. (15) 

For later use it is convenient to introduce instead of the initial angle ip a new 
variable 

u = VQsin- (16) 

which is suitable for the description of intermediate regions |(7J) where the usual 
diffraction coefficient description can not be applied. To investigate the multiple 
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scattering in such regions, we consider below the hmiting case Q oo with u 
kept fixed which corresponds to small initial angles of the order of 

(17) 



For reflected waves we shall distinguish between small angle reflection where the 
variable ^ 

Un = VQsin^ (18) 

remains finite when Q ^ oo and large angle reflection otherwise. 

In terms of variables H16|l and (|18|l the grating equation H10|l reads 



Un = \/n + v? (19) 

Physical values of reflection angles should obey the inequality < (p'j < tt or 
< sin((9j^ < 1. From it follows that the physical branch is deflned by the 
condition 

.W. = 2^^[l-^)>0. (20) 
In the notations of Fig. El the grating equation takes the form used in |Zj 



kp — uJnb — a^/k^-~LL;f^ = 2Trn, (21) 

where 

p — asinO + bcosd — dcosLp, (22) 
b — a/ tana is the shift of the half-plane corners along the z-axis, and 

uJn — k cos = k cos{a + ip'j^) (23) 

is the projection of the flnal momentum on the z-axis. The positive (resp. 
negative) sign of y^k^ — ojf^ correspond to forward (resp. backward) scattering. 

In particular, the specular reflection corresponds to n = in the above 
equation and 

uj() — k cos(2a — 9) = k cos(a + (p). (24) 
The possible transmission modes correspond to motion in a straight tube 

*(^™"")(z,x) = e'^^-^sin^mx, (25) 

where the transmitted frequencies are 



/TTTOX 2 , /_ / m 
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Here 

Q' = t^^Qsma (27) 

TT 

and 1 < TO < Q' is an integer. 
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3 Exact results 

At large distances the total reflected field is the sum over the individual allowed 
propagating modes © 

-Q sin^ (¥)/2) <n<Q cos2 (i/3/2) 

with ip'^ determined by Eq. pU|) and the coefficients Rn in the sum are reflection 
coefficients which deflne the probability amplitudes of reflection to the given 
channel. 

Similarly, the total transmitted field is the sum over individual transmitted 
modes (^3)) 

= J2 Trae^-^'-sin-mx (29) 

l<m<Q sill a 

with from Eq. H2()|) and the coefficients are called the transmission coef- 
ficients. 

We restrict ourselves to the case of the Dirichlet boundary conditions on the 
demi-plans treated in 7 . From the results of this work one can write down 
explicit expressions for reflection and transmission coefficients in the form of 
infinite products. 

The refiection coefficients in Eq. (I28|l with n = and n 7^ are slightly 
different. For n = the reflection coefficient Rq corresponds to the specular 
reflection ifi^ — ip and we shall call it the elastic reflection coefficient. From [7] 
one gets 



For non-zero n [7] 



The transmission coefficients in the expansion ()29|l are given by the following 
expressions 



(32) 



a^UJm{i^7n - " Wq) L ' J K{uJm) ' 

In these formulas cjq is given by 1)24(1 . w„ are deflned by H21|l and 

Lo' = k cos 9 — k cos(q; — (p) . (33) 
The function K{lo) is the ratio of two functions 

K{.) = f44e-, (34) 
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with 



.a,, TT 1 , /n • M 

X = —i — [{a — —) ^ — m(2 smajj, 



TT 2 tan a 

The denominator f+{Lo) in H34|l is given by the convergent infinite product 



(35) 



ui^) - n 



ka\ .oja 
■nn ) Tin 



(36) 



The numerator in (|34ll is the product of two functions 

g+{uj)^Gi(u;)G2[Lo), (37) 
where (71,2 (f^) are represented as the following convergent infinite products 



oo 



^(kp—ujb+iuja) /(27rn)+2(7r/2— a) 



n=l 
-1 



G2{uj) = Yl (A„ +z*„)e 



{kp—ub—iija) /(27rn) — i(7r/2 — a) 



and 



A„ = 



1 _ JieS] - ( Y 



' sin^ a 



27rn sin a 



27rn / 



(38) 
(39) 

(40) 
(41) 



The functions f+{oj) and g+{Lo) have no singularities or zeros in the upper half 
plane 

Im uj > Im Wo (42) 

where one assumes that the momentum k has a small positive imaginary part. 

These functions appear in the following Wiener-Hopf type factorization prob- 
lem 



sin[a\/^^ ~ ^'^] 



(43) 



and 
a2 + 62 



(a; — <jj^{u — LL)')g^{uj)g_{u)) = cos[a\/ k^ — u;^] — cos[fcp — ujb]. (44) 



Functions f-{oj) and g-{oj) have no singularities in the lower half plane 

Im Lo < Im Lu' . (45) 
The explicit form of these functions is 



/-(c^)=n 



, ka \ Lua 
1- — +i— 



Tin J nn 



-lauj / Tin 



(46) 
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and 

where 



oo 

n=l 

-1 

G2(w) = Yl (A„ -i*„)e('^'''-"''+*'^°)/(2^")+*('^/2-a)_ 



n— — oo 



The real zeros of the function g+{Lo) 



(47) 
(48) 
(49) 

(50) 



correspond to pure imaginary values of A„ and coincide with the reflected fre- 
quencies l|23f) . Their explicit form is 



k cos a cos If 



27rn 
~kd 



. 2 47rn / 7rn 
KsmaWsm 1^9 H — —- cos 1^9 — — 
kd V kd 



(51) 



In accordance with Eq. H20|l the positive branch of the square root has to be 
chosen. LOn with positive (resp. negative) n are real zeros of Gi{lj) (resp. of 

G2(W)). 

For the later use we need also the zeros lj* of the function g-{u}). They are 
given by Eq. (|51|l with the different sign of the square root 



k cos a cos {p — 

kd 



k sin ( 



47rn / ■nn\ 
— (cos^--). (52) 



aysm Lp - 



The transmitted frequencies uim are zeros of the function f^(uj). They are given 
by The zeros of /+(w) are just ~uj,n- 

The expressions (|30|) . (I31II . and (|32() for the reflection and transmission co- 
efficients include infinite products and are quite cumbersome. Simpler formulas 
can be obtained for the modulus of these coefficients which is sufficient for many 
purposes. Using repeatedly the defining relations (|43|l and H44f) one can demon- 
strate that the values of |i?„P and |T„p are expressed by the following finite 
products 



\R 



01 = 



n 



^0 - ^„ 



n 



LJo+UJ„ 



\Rn\ 



4 sin^ a sin^ 



(53) 



i^n - 




n 


- l^ni TT 






n 





f\\2 



(54) 
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4 sin tp 



and 



n 



'M2 



-Wo) 

n 



n 



/\\2 



\K{co')\ 



(55) 



\K{lo')\ 



n 



n 



OJ — OJ', U! 



(56) 



In Eds. (|53|I - H56|I the products are taken over finite number of real reflected (|51|l 
and transmitted H26(l frequencies. 

As w' — Wo = 2 sin a sin and ujq — uj' , Eqs. I|53(l and (|54|l can be rewritten 
in a form which is vahd for both ri = and n 7^ 



Wo 



n 



n: 



\K{co')\' 



(57) 



If 71 7^ the first product includes the term with n' = 0. 

The scattering amplitudes fulfill the current conservation. In the configura- 
tion of Fig. 121 it is convenient to consider the current through the surface of a 
rectangle enclosing the scattering plane whose two sides are parallel to it. Due 
to the periodicity of the demi-planes this surface is reduced to the one indicated 
by dotted line in Fig. [21 and the current through boundaries perpendicular to 
the scattering plane can be ignored. In this case the total current through the 
parallel parts of this surface is zero and a simple calculation gives the following 
relation between the reflected and transmission coefficients 



k sin ip — k sin p'^ \ Rn \ 



^sina I 

I 



fc2 



Trmy 
a / 



(58) 



All expressions in this Section are exact and suitable for numerical calculations. 
But for theoretical purposes they are practically intractable, especially in the 
semiclassical limit Q ^ 00, because the number of factors in Eqs. (|53|I - H55|I 
increases with Q and each factor depends non-trivially on Q. 

In the next Sections we investigate these expressions when Q 00. In this 
limit non-trivial results appear when the incident wave forms a small angle ip 
with the scattering plane (as in (|17|l ) and we focus our attention on this region. 



4 Elastic reflection coefficient 

Let us consider first the behavior of the elastic reflection coefiicient Rq in 
Eq. H30fl in the limit ip —> 0. Taking into account only the linear in p terms 
direct calculations give 

i?o = -(l + /?</?) (59) 
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and 



/3 = *E^" 



2 sin a 1 
V(l ~ Cn) ~ COS a (sin^ a + sin a) + z cos a 

1 



(sin a — sin a) — z cos a 

TT, 1 



-2z(3' sina[(c 



2"'' tan a 



— ln(2sinQ!)], 



where ^„ = Q' jn and Q' = kajn = Qsina. 

Separating the real and imaginary parts of Eq. H()()(l one obtains 



/3 = Q'(Ci+iC2), 



(60) 



(61) 



where 



[Q'l 



Ci = —2 sing 



1 VO^^+Q'cosr 



[Q] 

E- 



and 



Co 



[Q'l 
cos a ^ ■ 

n=l 



sma 



— cos a 



sma 



[Q] 

E 

=[Q']+i 



1 



?i — Q' sin a ' 



K-y/ (Q — n)n + n cos a 



(62) 



E 



yjn^ - Q'- 



E 



\/n + Q 



E 



n2_Q'2sin^a ^ \/"'("- + Q' sin a 

ji=[Q'] + l ^ n=l ^ ^ ^ 



ji=[Q] + l 



y/n{n — Q' sin a) 



-2[(a-5)-^-ln(2sina) 
2 tana 



(63) 



Here and below [x] denotes the integer part of x such that x = [x] + {x} and 
the fractional part {a;} obeys the inequality < {x} < 1. 

In the computations below we assume that a 7^ 0, 7r/2, and tt. In the 
semiclassical limit fc — > c» all sums in the above expressions include many 
terms which can be estimated by usual methods. The dominant contribution 
corresponds to the change of the summation to the integration. One gets 



-2 sin a 



Q' 



dn 



Jo 



dn 



yQ^ — r? + Q' cos a Jo sin a \/ {Q — n)n + n cos a 



cos a 



Q' 



dn 



n + Q' sin a 



cos a 



dn 



n — Q' sin a 



(64) 
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and 

C2 = sin a 



-an — / — --— rdn 



Qi — (5'2 sin^ a Jo Vn{n + Q' sin a) 

^ (65) 



- r ri'^?- ^rfn-2[(a-;)-^-ln(2sina)] 
Jq ysmaj z tana 

The integrals are elementary and cancel each other, i.e. Ci = C2 = 0. One can 
prove that the first correction term (of the order of l/\/Q) to Ci appears from 
small n summation in the second term in Eq. I|62(l 

[Q] ^ 

Ci 2_, — h smooth terms 

sin a^/ {Q — n)n + n cos a 

^ hm > ^ - / ^ (66) 

sina^/Q V" Jo J 

The integral is subtracted because we know that all integrals over n are canceled 
by other terms. 

The last limit can be computed from the relation 



N \ 1 

= CU (67) 



Af-^oo \ ^ yfn 2 2 

\n— 1 ^ / 

where C(s) is the Riemann zeta function ~ —1.460354). 

Finally when Q — > 00 

ci - -^ah- (68) 

smavQ 2 

Similarly the dominant contribution to C2 comes from small-n summation in 
the second term of Eq. (|63|l 



00 



C2 ^ - sina V ZV^- V + smooth terms ^ =C(i). (69) 

Vrtin + y sma smavO 2 

From Utilfl one concludes that for large k and a 0, 7r/2, n 

/3=^(l_z)C(i). (70) 

In terms of variable (|16|l this result states that the two first terms of the expan- 
sion of the elastic reflection coefficient into the power of u are the following 

Ro = -l-2V2e-^^/\{^)u. (71) 

In Section |^ we demonstrate that these two terms can be obtained from the 
KirchhofF approximation developed in 
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5 Small angle reflection in the limit ^ 

When the incident angle it follows from Eqs. (|31|l and H32|l that the re- 
flection coefficients with n ^ and all transmission coefficients are proportional 

to ip. 

In the computation below it is convenient to rescale all frequencies by k i.e. 
to change 

To simplify the notation we shall from now use the same symbols for the rescaled 
quantities. 

When ip ^0 rescaled Eqs. (IHU, lE^J, and give 



cos«(l-|)-2sinJg(l-^), (73) 



cosa(l-|)+2sinJ^(l-g), (74) 



and 



y sm a 

with all square roots chosen positive. 

Using Eqs. H54|l and 1)56(1 one can compute the reflection coefficient modulus 
in limit fc cxd by writing the logarithm of each product as a sum over the 
corresponding frequencies. E.g. 

log |cj„-a;„'| = ^ log|a;„-w„'|. 

n' n' 

When k oo the sum can be substituted by the integral. Exactly as it was done 
in the precedent Section one can check that all integrals in the full product H54|l 
cancel and the dominant contributions come from regions with small factors. 

We are interested first in the behaviour of the reflection coefficients i?„ at 
fixed n ^ 0. In this case there are 3 regions with small factors. The first 
corresponds to u;„ — ujrn' — > 0, the second appears when a;„ — w*, 0, and the 
third includes cases when n' is close to n (see Fig. (Ol). In the first region one 
should expand the factors near Q sin^ a, i.e. 

n' = Q sin^ a + Sn, (76) 

where 

Sn = -{Q sin^ a} + q (77) 

with integer q and the similar expansion for to'. (As above {x} denotes the 
fractional part of a;.) 
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Figure 3: Dependence of different frequencies on n (0 < n < Q). Left: < a < 
7r/2, right: 7r/2 < a < it. w„ is represented by solid line, by dotted line, w„ 
by dashed line, and — cD„ by long-dashed line. The full small circles indicate the 
point with coordinates (Q sin^ a, cosa) where two different frequencies (with 
^p — Q) coincide. 

The required series when < a < 7r/2 are easily obtained from Eqs. (jZSJ- 



^ (fa) ^ Q (^(fa^3 

^cosQ! 4(52 sin^ a cos3 a \ Q 



and 



5m {5raY , /nfi-^^\3 



UJm\m=n'+5rn = \ COS a\ - — - \ — + O ( — )^ . (79) 

Qlcosal 2(32 sm^ q;| cos^ a| \ Q J 

For completeness we add 

Wn|n=n*+5n = COSa(l - 4 siu^ a) + 0{^). (80) 

At small n <^ Q 

/Tl TL I Ti Tl 

— +C'( — ), w* = cosa + 2sina. /— + C'( — ) (81) 

In the first region when n is small and < < 7r/2 one formally has 

Similarly in the second region the expansion of lj„ — cj*, gives 

*^^TT/^ o- -{Qsin^a} +q\ 
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The two products (|82|l and H83|l are identical and cancel each other in the ex- 
pression for the reflection coefficient. The same remains true when 7r/2 < a < ir. 
Therefore a dominant contribution comes only from the region of small n' where 



" 1 



71' n'^7i V n 

exp(-/ [fog(l+ /Il)-log(l- (84) 



n' V n' 



The last term is added because the total integral is canceled by other terms. 
When ^ cx) the integral equals + 0{l/N). 

Introducing the convergence factors one gets that this product tends to the 
following function 



exp 2v/7^(;^-^-^/]v/2) \[ 



1 - 



71 



e^'v^. (85) 



When — > oo using Eq. 167|l we obtain 



n ^71 



Combining all factors together one concludes that in the limit Q — s- cx) the first 
term of the expansion of the modulus of the reflection coefficient at small n ^ 
has the following asymptotics 

|i?„r^°g^'|r„P=4^.2|r„|2, (87) 
where u is defined in (|16fl and 

|r„P = ^. (88) 



In Fig. ^we present the results of numerical calculation of the first term of ex- 
pansion of reflection coefficient from Eq. H54|l together with asymptotic formula 
(|88|l for this quantity. The agreement becomes better at larger Q as it should 
be. 

Notice that (when a is not too close to 0, 7r/2, and tt) the reflection coef- 
ficient at small n do not depend on a. One can check that when Q — > cxd the 
transmission coefficients with small m are always negligible. 
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50 



Figure 4: The first term of the expansion of the reflection coefficient at small 
n > 1 into a series of u. Two upper curves are |r„p computed at Q = 1000 and 



Q = 8000 with a 
this quantity. 



7r/3. The lower curve is the asymptotic expression for 



6 Reflection and transmission coeflicients at 
large n and cp ^ 

One can prove that for large n the only important contribution comes from 
n close to n* = Qsin^a where two different frequencies coincide (see Fig. (S)). 
When < a < 7r/2 in this region only the transmission is noticeable and when 
tt/2 < a < TT only the reflection is important. 

Let us consider first the case < a < 7t/2. The expression (I55II for the 
transmission coefficient modulus |T„ip can conveniently be rewritten in the 
following form explicitly separating small factors 



|T„ 



. 2 

sm ip 



n 



UJr, 



n 



n 



u' + LOri 



n 



UJ — LOrn.i 



UJr, 



(89) 



In these formulas n' and m! are non zero and we assume that m = Q siv? a + 5m 
where 5m is of the order of \/Q. 

From Eqs. f7H| . (|7^ . and (P?T|l formal expansions for the above products can 
be obtained. The product over n' contains terms with small n' — n and large 
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n' = Qsin^ a + 5n' . Putting 5n' = 5m + q and taking into account dominant 
terms one gets 



7i' n>l 




(90) 



The first product comes from small n and the second one is due to the contri- 
butions with large n' = Q sin^ a + 5m + q. 

Computing other products in the similar way one obtains that the first term 
of the expansion of the transmission coefficient into a series of u has the form 

i7^™r-" 77^5(«/), (91) 

Q sm 2a 

where the scaled variable m/ is related to m as follows 

m ^ Q sin^ a 



Q sin 2a ' 



(92) 



and the function g{x) is 



, 2 

X 



/n 

The renormalization factor e^''*^^^^ has been computed exactly as it was done 
above for the reflection coefficient at small n. 

When 7r/2 < a < tt one can similarly show that transmission coefScients at 
small ip are negligible and the reffection coefficients with n close to n = Q sin^ a 
differ only by the factor 1/2 from the transmission ones with the same value of 
n and a = tt — a 

K(a)r- ^|r„(7r-a)p. (94) 

For clarity we explicitly show the dependence of a. The factor 1/2 with respect 
to Eq. H91|l is related with the different normalization of real transmitted modes 
and complex reflected ones ©■ 
The formulas above are valid when m/ is fixed and Q ^ oo. The meaning of 
the variable uj (|92|) is most easily seen when 7r/2 < a < tt. The scattering demi- 
plans can be considered as a system of mirrors and the initial ray with = 
after the specular reflection in these (demi) mirrors will form the angle 2n — 2a 
with the scattering plane. For the reflection close to this optical boundary (cf. 
©) it is naturally to introduce a new variable ?/;„ which measures the deviation 
of the reflected angle from the angle 27r — 2a 

Lp[^ = 27r - 2a - (95) 
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Figure 5: The first term of the expansion of the transmission coefficient for 
Q — 1000 and a = 7r/3 into a series of u (dots). The solid Une is asymptotic 
formula H93|l 

At large Q the variable in (|92() (with m — n) is related with ^„ exactly in 
the same way as the variable u is related to the incident angle (see (jl6|l ^ 

Uf = ^/Qsu\^. (96) 

We write Uf without the index n in order to stress that for large angle scat- 
tering the discreetness of reflected and transmitted waves plays no role and the 
asymptotic formulas H91|) . (|94|l can be considered as functions of a continuous 
variable Uf. 

In Fig. |3l we present the results of numerical calculations of the function 
g{uf) calculated from Eq. H55(l together with the asymptotic formula H91|l . The 
agreement is quite good. The reflection coefficient for 7r/2 < a < tt is also well 
described by the asymptotic relation (|94|l . 

7 Reflection and transmission coeflicients at fi- 
nite u 

Eqs. (|87|l . (|91|l and (|94|l represent the first terms of the expansion of the trans- 
mission and reflection coefficients in the power of u = -y/Q sin(</j/2). 

The purpose of this Section is to calculate these coefficients at finite values 
of u. The calculations are performed exactly as in the precedent Sections. First 
we note that the frequencies LOn and cj* depend on u only in the combination 
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Figure 6: The reflection coefRcients with n = —3, . . . , 3 with Q = 8000 and 
a = tt/3 (sohd lines) together with the asymptotic formula H97|) (dots) for the 
same values of n. 



n{u) = n + u^ . Therefore precedent formulas remain valid when n is substituted 
by n{u) = n + u^ and the dominant contributions come from the same terms as 
above. We omit the details and present only the final answers. 

The reflection coefficients at all small n (for negative n > — [u^], n — 0, and 
positive n) are given by 



|i?„(u)|' = 



2u 



(97) 



and the function G{x) has the following form 



G{x) 



n 



_9_£_ 



n 

n'>0 



(98) 



The terms with the same value of n' should be grouped together and then the 
total product converges. 

In Fig.Elthe results of numerical calculations for |i?„(u)p with n = —3, ... ,3, 
Q = 8000 and a = 7t/3 are compared with the asymptotic formula (|97|l and the 
excellent agreement is found. 

In the discussion of the scattering at large angles corresponding to n close to 
Q sin^ a it is convenient to introduce reduced transmission t{uf,u) and reflection 
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r{uf,u) coefficients in the following way 



and 



T — 



Q sin 2a 
2 



r{uf,u), 



(99) 



(100) 



/Q sin 2a 

where the variable u defined by Eq. Hlt)|) fixes the initial angle and the variable 
u/ is related with the scattering angle by Eq. (|92|l . 

As above one can demonstrate that when < a < tt/2 only the transmission 
coefficients are important and when 7r/2 < a < tt the reflection coefficients 
dominate and asymptotically in corresponding regions 



|r(u/,u)p = -|t(M/,U,)p ^g{uf,u) 



(101) 



where 



y2g2(«+«y)C(i) 



n 

n'>[«2] + i 



n 



n 

n'>0 



(102) 



As above, the terms with the same n' should be combined together for conver- 
gence. 

The reduced coefficients t{uf,u) and r{uf,u) introduced in Eqs. (pl^ . (|100|l 
can be considered as the reflection and transmission coefficients for the scat- 
tering in the (continuous) interval of 'angles' Uf as it follows from the current 
conservation H58|) rewritten in terms of new variables 



/oo 
dufg{uf,u). 
-oo 



(103) 



n> — \u 



Here i?„(u) mean the reflection coefficients at small n. 

When u ^ oo the precedent equations can considerably be simplified by 
noting that in this limit the products in these equations are large when is 
close to — li. Putting u/ = —u + 5 and taking into account the dominant terms 
one gets 



g{uf,u) 



sin(27r(uy + u)u) 



(104) 



27r(u + Uf) 

Practically this approximation works well even with u > 1. The semiclassical 
derivation of this expression is performed in Section |^ 
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Figure 7: Reduced transmission coefficients with a = 7r/3 (left) and the reduced 
reflection coefficients with a = 27r/3 (right). In both figures u = 2.1. Dashed 
line: Q = 1000, dotted fine: Q = 4000, sofid fine: Q = 16000. The thick solid 
line are the asymptotic formulas (|102|l and H101() . The thick dotted line is the 
approximation (|104ll . 

In Fig. we present numerically computed transmission coefficients with 
u = 2.1 and a — vr/S and reflection coefficients with a = 27r/3 for Q — 1000, 
Q = 4000, and Q = 16000 together with the asymptotic formulas H1U2|) and 
P01(l . With increasing Q the agreement becomes better and better. The ap- 
proximation H104() is hardly distinguished from the exact asymptotic formula 

8 Scattering when a = 7t/2 

The formulas of the precedent Sections are not valid when a = 7r/2, i.e. when the 
demi-plans are perpendicular to the scattering plane. Nevertheless, asymptotic 
expressions for large Q can also be obtained in this case but the results depend 
on the fractional part of Q and, strictly speaking, the semiclassical limit Q — > oo 
does not exist. The main reason for such behaviour is related with the existence 
of a trivial solution for some particular values of Q. It is easily seen that when 
a = tt/2 the following function (in the notation of Fig. [21) 

^'(0) (a;, y) = sin(fc2; cos ^)e*'''^ '^'^ (105) 

is an exact solution of our problem (i.e. it vanishes at the demi-plans) provided 

kdcosip = ttI, (106) 

with integer /. In the notations (|13() and H16() this relation takes the form 

{g-2u2} = 0, (107) 
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where < {x} < 1 is the fractional part of x. 

Irrespectively how large is the dimensionless momentum Q, close to points 
when Eq. H1U7|I is fulfilled, the appearance and disappearance of the exact so- 
lution H1U5|I strongly perturb the reflection and transmission coefficients. 

Though the limit Q oo does not exist, a special limiting case when {Q} 
fixed and [Q] — s- oo can be computed from the approach developed in the pre- 
ceding Sections. We omit the details and present only the final formulas. 

Let us define a function W{t) as follows. 



W{t) = e2(2-v^)tC(i) -Q ^-'^'-^''7^ 




(108) 



We shall need this function at special values of the arguments t — y^n + {u^}, 
t — \J n + {Q — u^} and t = — -^-^/n + {Q} with integer n. The prime in the 
second product in Eq. H109() means that when t = ■y/n + ju^} the term with 
n' — n IS. omitted in the first factor, when t = \/ n + {Q — m^} the term with 
n' = n is omitted in the second factor and when t = —^=^n + {Q} the term 
with n' = n is absent in the third factor. 
The refiection coefficients at small n 

\Rn{u)\^ = f W{u^)W{u), (109) 

where m„ — \Jn + v?- determines the allowed value of small reflection angle (see 
(O and IHl). 

The refiection coefficients at large n — \Q\ — g is given by the same expression 
Hl()9|l but Un is substituted by Uq 

where 

Uq = ^q + {Q}-y? (Ill) 

has the meaning of the allowed value of small deviation of reflected angle from 
its maximum possible value. It means that if one writes Lp'^ = tt — 5ip'^ then for 
n^[Q] — q with large [Q] and fixed q 

sin^^Mzi. (112) 
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Figure 8: Reflection coefficients witfi n ~ —3, ... ,3 witli a = tt/2 and Q = 
8000.7 (solid lines) together with the asymptotic formula (|109|l (dots) for the 
same values of n and {Q}- 



The reflection coefficient Rq with integer g = 0, 1, . . . , exists for a finite interval 
of u such that < w < y/q + {Q}- 

The transmission coefficients for m = [Q] —p with p = 0, 1, . . . , are given by 
the same expression (|109|l but with the substitution u„ Up 

\TM? - i^^]^ W{up)W{u), (113) 



where 



If for a — -K /2 one determines the angle of transmission 0m from the natural 
relation cos ipm — m/Q then for m — [Q] — p Up is related with (5(/)p = tt — 
4'm=[Q]-p ill same way as above 

Q,-J_MzL, (115) 



Ay Figs. 181 101 we present the the results of numerical calculations of the 
reflection and transmission coefficients together with the above asymptotic for- 
mulas. The agreement is very good. 
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Figure 9: Elastic reflection coefficient with a = tt/2 and [Q] = 8000 for 
different values of fractional part of Q (solid lines). From top to bottom 
{Q} = .001,0.201,0,401,0.601,0.801. For clarity each curve is lowered with 
respect to the precedent by 2 units. Dots represent the asymptotic formula 

in 




Figure 10: Reflection coefficients with a = 7r/2 and Q = 8000.7 for the largest 
values of n = [Q] — q with q = 0, . . . ,3 (solid lines). Dots represent the asymp- 
totic formula H110|) . 
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Figure 11: Transmission coefficients with a — ^[2 and Q — 8000.7 for the 
largest values of m = {Q}—p with p = Q and p = 3 (solid lines). Dots represent 
the asymptotic formula \\V6\ . 

9 KirchhofF approximation 

Semiclassical limit of different cases of multiple diffraction near the optical 
boundaries have been considered in Ref. [Sj. In particular in this paper the 
contribution to the trace formula from diffractive orbits close to n-fold repeti- 
tion of a primitive diffractive orbit has been calculated. 

In the Kirchhoff approximation this contribution is given by the diagram of 
Fig. El Each line in this figure corresponds to the free Green function 

Go(f,f) = — ==^, (116) 

where I = \x — x'\ is the distance between two points af, x' and k = ^/E is 
the momentum. Each circle at Fig. El describes the convolution of two Green 
functions and in the Kirchhoff approximation it gives the factor —2ik. The role 
of the obstacles (corners) consists in the restriction of the integration over yi to 
the half line. 

In Ref. in] it was demonstrated that the contribution to the trace formula 
from such trajectories has the form 

pW/)(£;) = _ J.^A„e'''*' + C.C., (117) 
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Figure 12: Schematic representation of multiple diffraction near optical bound- 
ary of n equally spaced obstacles. 



where An is given by the following n-fold integral 



A-n — 

with 



i/2 



dsi 



ds 



2 • ■ 



dSn 



dso 



$(S ) = (Sl - 32^ + (S2 - Ss)^ + . . . + {Sn-1 " S^f + (s„ 



dSn 



(118) 



(119) 



The presence a discrete symmetry in this quadratic form permits the analytical 
computation of this integral and |S] 



-J n— 1 



1 



When n oo the sum over q can be substituted by the integral and 



lim A,- 

n — >oo 



1 



dq 



71- Jo y/q{n - q) 



1 



(120) 



(121) 



and as it was noted in Ref. 6. the contribution 1117(1 in this limit coincides with 
the contribution from the boundary trajectory reflected from a straight mirror 
with the Dirichlet boundary condition. 

Let us consider this point in details. Assume that after each reflection with 
a mirror a trajectory gets a reflection coefficient R. Then its contribution to 
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the trace formula computed in the usual manner is 



ifcL — 3x1/4 rL ji 
^{boundary) ^^-^ ^ -2i?^ / " / J^yg^'^ky'' {1/1 + 1/ {L-l}) 



IGnVk Jo y/l{L - /) Jo 

Here L is the length of the trajectory and the factor 2 takes into account that 
each trajectory can be passed in two directions. For simplicity we do not consider 
the symmetry factor and consider the trajectory as being primitive. 

Comparing this result with Eqs. pi7|l and (|121|l one concludes that these 
equations can be interpreted as the reflection from a mirror with the effective 
reflection coefflcient i? = — 1 as from a boundary with the Dirichlct boundary 
conditions. 

To compute the behaviour of the reflection coefflcient for large but flnite n 
it is convenient to use Eq. (|67|l from which it follows that when n oo 

A„^l + ^^. (123) 

One can incorporate this result into the above picture of the reflection from 
a straight mirror by assuming that the reflection coefficient R depends on the 
reflection angle ip (or transverse momenta py « kip) 

R{e) = -{1 + P^). (124) 

As in the saddle point approximation ip w "iy/L, the only modification of the 
above calculation is the following integral 







Because L — In one concludes from (|123|l that 

V TT z 



(126) 



which up to notations agrees with Eq. (|71|l obtained by the direct expansion of 
the exact solution. 

The above considerations demonstrate that small-angle singular refiection 
from a periodic set of demi-plans can be interpreted as much simpler process of 
the specular refiection from a straight mirror whose reflection coefficient at small 
angles has behavior as in Eqs. (|124|l and (|126|l . Of course, at small distances 
reflection flelds for two processes are very different but at large distances they 
are equivalent. 

Similar (but simpler) considerations permit also to understand the approxi- 
mation (|104|l for the reflection (and transmission) at large angles. 
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reflected rays 




Figure 13: Schematic representation of the scattering at large angles. Black 
circles are the corners of the scattering demi-plans. The dotted line represents 
the scattering plane. 

At Fig. ^] we represent schematically the configuration important for the 
large-angle scattering. The rays reflect from small parts of the scattering demi- 
plans restricted by the indicated corners. The width of each effective mirror 
is 

A = dsinv?, (127) 

where, as above, d is the distance between singular corners along the scattering 
plane and ip is the scattering angle. After unfolding the amplitude of such 
reflection can be calculated in the Kirchhoff approximation as the transmission 
through a slit of the same width by the usual formula (an additional minus sign 
is due to the reflection from a mirror) 

D{e) = 2ik / e-'''y ^""dy = sin -fcA sin e , (128) 

J -A/2 sme V2 / 

where e is the deviation of the scattering angle from the direction of the specular 
reflection. 

The total reflected field is the sum over all diffracted fields 

oo ikR — 37r't/4 

^^^^f\x,y)^ V e'^'^™™^^Z?(V)^7r^f^, (129) 

where R„i — ^^(x^^^Tndp^r^ is the distance between the m^^ diffractive cen- 
ter and the point of observation with coordinates {x,y). Using the Poisson 
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summation formula one gets 

°° /-oo -37ri/4 



(130) 



where 

S{m, n) = —2Trmn + k{dm cos tp + Rm)- (131) 

When k ^ oo one can use the saddle point method. The saddle point is obtained 
from the condition dS{m, n)/dm = which gives the grating equation pU|) . The 
computation of the second derivative and the resulting integral leads to 

\]/(''e/)^2; y) — \^ Dj^n) ^ik(x cos ip'^+y sin ip'„) (132) 
n 

where the reflected angle is defined by Eq. (fTU)! . 

Using Eq. (|128|l for the diffraction coefficient of the reflection on a small 
mirror and taking into account that in the important region (p'^^ k.Iix ~1a one 
obtains the following expression for reflection coefficient _R„ 

2 1 

i?„ = ---— — - — : sin(-fc(isini^sine„). (133) 

kd sm la sm e„ 2 

As (/5 ~ lujyfQ and e„ = —2{uf + u)/^/Q with Q = kd/n, the last expression 
equals 

i?„ = r sin2^.(., + u) 1 ^ ^^3^^ 

whose modulus coincides with Eq. (| 104(1 . 

10 Summary 

The exact transmission and reflection coefficients for the scattering on infinite 
number of parallel demi-plans obtained in Ref. [7] are analyzed in semiclassical 
limit of large momentum. 

The most interesting (and difficult) case of small incident angle is consid- 
ered. More precisely, the incident angle ip is chosen in such a manner that 
approximately (p = 2u/^/Q where Q is the dimensionless momentum. The limit 
considered corresponds to Q — > oo with fixed u. 

It is demonstrated that at small final angles (of the order of Un/VQ) the 
transmission is always negligible. The modulus of reflection coefficients in this 
case are independent of Q and the angle a of inclination of the demi-plans and 
are given by Eq. (|97() . The reflection coefficients decay quickly with n (i.e. with 
increasing of reflection angle) and J2n \/"-|^nP converges. The largest reflection 
coefficient corresponds to the smallest possible angle of reflection i.e. n = — [w^]. 
For very small incident angle the elastic scattering corresponded to the specular 
reflection dominates. 
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The large angle scattering is noticeable only close to the specular reflection 
from the demi-plans. When < a < 7r/2 only the transmission is large and the 
large-angle reflection can be neglected. The transmission coefficients have the 
form t{uf,u)/{\/Qsin2a) where \f{u,Uf)\^ = 2g{uf,u) and g{uf,u) is given 
by Eq. (|l()2|l . Here u/ is the deviation from the angle of mirror reflection 
magnified by the factor exactly as u is related with the incident angle. 

For tt/2 < a < TT the transmission is small and the reflection coefficients have 
similar asymptotics. 

The exceptional case of demi-plans perpendicular to the scattering plane is 
characterized by the dependence of the fractional part of the momentum. In 
the limit {Q} fixed and [Q] — s- oo the reflection and transmission coefficients are 
independent on [Q] and are given by Eqs. H 109(1 - H113() . 

The first two terms of expansion of the exact elastic reflection coefficient into 
powers of small incident angle can also be obtained from the results of Ref. [H] 
where the Kirchhoff approximation for multiple scattering was developed. The 
large angle scattering is well described by the usual Kirchhoff approximation 
when the initial parameter u > 1. 
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